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framework developed in earlier work [l|, ^, which showed how to associate entropy currents 
with horizons in the dual geometry. The entropy current defined by the event horizon in 
the dual bulk geometry is calculated. It is also shown that to second order in the gradient 
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result obtained on the basis of a purely hydrodynamic analysis |Q. 
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1. Introduction 

For static black holes in Einstein gravity there is a well established notion of entropy 
expressed in terms of the event horizon area by the Beckenstein formula Q]. One would 
expect that this notion should, in some form, apply also to black holes which are in an 
appropriate sense close to being static |^. However, the event horizon boundary of a 
black hole is a global, teleological concept ~ its location can only be determined once 
the complete spacetime is known, since one must examine the ultimate fate of all signals 



originating from a point before ruling whether or not that point is part of the black hole 
Q. In this way the evolution of an event horizon can appear to be acausal. The area of the 
event horizon is similarly teleological. For example the area increase of an event horizon is 
not directly driven by infalling matter or energy; the actual effect of an influx through the 
event horizon is to decrease its rate of expansion [^]. This becomes a puzzle as soon as one 
wishes to interpret the area of the event horizon as a measure of black hole entropy, since 
the apparently acausal expansion of event horizons would then seem to imply a similarly 
acausal evolution of entropy. An acute expression of this comes arises in the context of the 
AdS/CFT correspondence, which maps the entropy of an AdS black hole to the entropy of 
a fluid in the dual quantum field theory on the conformal boundary. 

It is therefore an open question whether a naive reading of the Beckenstein formula accounts 
for black hole entropy for dynamical black holes. The simplest option is to assume that 
the the entropy is still proportional to the area, but perhaps not to the area of the event 
horizon itself, but of a hypersurface which asymptotes to it at late times [^, ^, |l|, ^. One 
such example is the apparent horizon, understood as the boundary of a region containing 
trapped surfaces. For static black holes the apparent and event horizons coincide, but the 
location of an apparent horizon on a given slicing of spacetime can be determined locally 
in time. Since the apparent horizon evolves in a causal way the associated entropy also has 
this property. Such a notion, while not free of conceptual problems, has attracted a lot of 
attention Q. 

In the context of the AdS/CFT correspondence the long wavelength distorsions of black hole 
(or black brane) horizons map to hydrodynamics of the dual quantum field theory plasma 
M. The notion of near-equilibrium entropy then maps to the hydrodynamic generalization 



of entropy - the entropy current |1C, 11 1. Such a current is constructed phenomenologically 
in the gradient expansion by requiring that in equilibrium it reproduces thermodynamic 
entropy and that its divergence evaluated on solutions of the equations of hydrodynamics 
is non-negative. A detailed analysis of the consequences of this generalized second law 
of thermodynamics on the form of the entropy current in [12, ^] showed that already at 



second order in gradients there is an ambiguity inherent in such a definition. 

It is very natural to suspect that the bulk counterpart of this ambiguity is related to the 
choice of horizon assumed to be the carrier. In the case of hydrodynamics with no conserved 
charges beyond the energy momentum tensor it was shown in [|^, y] that this is indeed the 
case. This was done by providing an explicit formula for the entropy current associated 



with a given hypersurface satisfying the area theorem^ . This formula reproduces the result 
of [llOi when evaluated on the event horizon. 



From the perspective of the phenonienological definition of the hydrodynaniic entropy cur- 
rent none of these hypersurfaces and none of the available bulk-boundary maps is favored 
over any other. However, causality of the boundary field theory seems to favor the entropy 
current dual to the apparent horizon - provided that it is free of the ambiguities related 
to foliation dependence and that the bulk-boundary map in use is causal. These issues 
discussed at length in Q. The present article applies the same methods to a more com- 
plicated, but practically important case of hydrodynamics with a conserved current. The 
dual holographic description of this system was established in |15,[l^ (see also |17, 18|). In 



this paper the solution appearing in [19| is used, which is somewhat more general in that 



it allows for a weakly curved boundary geometry and is presented in a gauge convenient 
for the task at hand. 

The organization of the paper is the following. Section |2| briefly reviews the geometry dual 
to conformal fluid dynamics with a conserved charge current obtained in [19|. Section ^ 



describes the calculation of the relevant apparent horizon in this geometry up to second 
order in gradients. Section |^ applies the formula introduced in [|l| to the case at hand. The 
entropy currents defined by the event horizon and by the apparent horizon are found and 
compared. A general discussion of the results and possible future directions of research is 
provided in Section 0. Appendix ^ provides a very brief account of Weyl invariance, and 
appendix y lists Weyl invariant tensors needed in the main text. Appendix ^ provides some 
details of the second order hydrodynamics which follow from the solution of ||l^ . Finally, 
appendix y| contains explicit results for the expansions needed for the computation of the 
location of the apparent horizon. 



2. The geometry of fluid-gravity duality with a conserved charge 

To describe hydrodynamics with a conserved current a Maxwell field in the bulk is required 
|15, 16, 19|. The action of the five-dimensional Einstein-Maxwell theory under consideration 



^A formula of this type is also discussed in |l4[] . 



reads^ 

where Gj\f is the 5-dimensional Newton's constant and the cosniological constant is denoted 
by 12/L^. For this theory to be a consistent truncation of type IIB supergravity[20, 21, 22 
the Chern-Simons couphng k has to assume the value l/2-v/3. 



The equations of motion derived from ( ^.iD support a 4-parameter family of exact, static 



black hole solutions |15, 16, 19] with planar horizons obtained by boosting and dilating 
the AdS-Reissner-Nordstrom black brane solution [^] . The constant dilation parameter is 
denoted by b and the boost parameter u^ is a 4-component velocity vector in the "boundary 
directions"^. The solution can be expressed in the form 

ds^ = r^ (P^^ - 2Bu^u^) dx^'dx'' - 2u^dx^'(lr, (2.2) 

where^ 



^=2V^ b^ 



Jz(l + ^^^^)+^) (2-3) 



and 



P^lv = h^y + u^Uy (2.4) 



is the projector operator onto the space transverse to u^, The vector potential takes the 
form: 

A = ^n,dx^ (2.5) 

The geometry ( |2.2D may be regarded as a stack of constant-r 4-dimensional planes, starting 
from the conformal boundary at r = cxd (which is 4-dimensional Minkowski spacetime) 
down to the curvature singularity at r = 0. The latter is shielded by the event horizon at 
r = 1/6. The parameter b appearing in ( p.3| ) is related to the Hawking temperature T of 
the event horizon by 

T=^^{2-q%^). (2.6) 

The lines of constant x^ in (|2.2| ) are ingoing null geodesies, for large r propagating in the 
direction set by u^, and the radial coordinate r parametrizes them in an affine way |1C]. 



^Latin indices run over all spacetime dimensions, i.e. from to 4, while Greek indices run over the 
"boundary directions", i.e. from to 3. 

^The 4-velocity vector normalized so that m^u*^ — —1 in the sense of the metric on the conformal 
boundary of the locally asymptotically AdS spacetime (p.2|), whose components are denoted by h^v 

^The notation is chosen so that in uncharged limit (g — >■ 0) B is equal to B{br) as defined in ||23|| . 



If b, u^ and /i^jy are allowed to vary slowly compared to the scale set by b, the metric (^ 
should be an approximate solution of nonlinear Einstein's equations with corrections orga- 
nized in an expansion in the number of gradients in the "boundary directions" parametrized 
by x^. This gradient-corrected metric can be written in the form 

ds^ = (g^^ _ 2n^V^) dx'^dx'' - 2u^dx'' (dr + rA^dx") , (2.7) 



with the condition u^Q^u = completely fixing the gauge freedom [23|. The field A 
appearing in ( |2.7D is the "Weyl connection" defined in appendix |A[ It is clear that lines 
of constant x are geodesies (affinely parametrized by r), as in the case of the static metric 
(|2l 



The form of the metric ( |2.7| ) is restricted by the Weyl-invariance^ of the bulk theory |^3[ : 
the functions V^ are of unit Weyl weight and Qfj_u are Weyl invariant. They can be expressed 
as linear combinations of independent Lorentz vectors and tensors built out of b, u^, h^y 
and their derivatives, order by order in the gradient expansion. At first order in gradients 
one has a Weyl-invariant vector 

If, = e^^xpu'V^P (2.8) 

and pseudovector 

Vo^ = q-^P;V,q (2.9) 

and a second order symmetric tensor^ 

cr^^u = -jT^ipUu) (2.10) 

of Weyl weight —1. The independent objects appearing at second order in gradients are 
listed in appendix |^. 

As shown in fl^, the resulting solution, up to second order, takes the form 

6 5 

V^ = r^Bu^ + rFi/^ + br^FoV^^ +r^Y. ^^^^^/^ + ^ Y. ^^^^M' 

6 11 

Q^,u = r^P^^u + 2br^F2(J^u + r^Y. ^iSiP^^u + J2 ^iTt^u- (2.11) 

i=l i=l 

The 31 coefficient functions Fq,Fi^F2, Ki,Li {i = 1,...,6), Wi {i = 1,...,5), and Hi (i = 
1, ..., 11) all depend on the Weyl invariant variables b'^q and br. 



A brief definition of Weyl invariance appears in appendix |A|. 
^Symmetrization is defined as A(^i,) := A^u + A^^. 



The Maxwell gauge field ^ is a vector field of Weyl weight zero. It will be taken in the 
gauge Ar = and its form is [|l9l 



6 5 \ 

+ Yolf, + YoVo^ + r Y, NAu^ + J^ Y,V^^ dx^ . (2.12) 

1=1 i=i J 



Again, the coefficient functions Yo,...,Yq and Yq and Ni {i = 1,...,6) depend on the 
variables br and h^q. 



All the coefficient functions appearing in ( p.ll ), ( 2.12| ) have been determined by solving 
the equations of motion^ and are given explicitly in [jl9[ . 

The metric given above is a solution of Einstein equations with negative cosmological con- 
stant up to second order in gradients, provided that h and u^ satisfy the equations of 
hydrodynamics, i.e. the equations of covariant conservation of the energy-momentum ten- 
sor and charge current obtained from ( p. Ill ) by holographic renormalization [^]. Explicit 
formulae |19] can be found in appendix (^. 



3. The Weyl-invariant apparent horizon 

This section is devoted to locating the Weyl-invariant apparent horizon for the spacetimes 
defined by the metric (^]^) following the approach of |j^, ^]. The time-evolved apparent 
horizon denoted by A is presented as the level set of a scalar function S{r, x). This function 
is required to be Weyl-invariant and can be written in the gradient expansion in terms of 
all the independent scalars available up to some order. Since the event horizon at order 
zero is at r = 1/6 and at that order A should coincide with the event horizon, one has 

S{r,x) = b{x)r — g{x) , (3.1) 

where g{x) is a Weyl-invariant function expanded in gradients of b, q: 

9{x) = go{x) + gi{x) + g2{x) + . . . (3.2) 

Here g^ denotes a linear combination of all Weyl-invariant scalars at order k in the gradient 
expansion. Thus go{x) is a function of Weyl weight zero, i.e. it depends only on the 
combination b^q. There are no Weyl-invariant scalars at order 1, and 6 at order 2, so one 



''This, as well as most other computations in this paper, was done with the help of Mathematica. 



expects to find 

go{x) = X{b^q) , 

9i{x) = 0, 
6 
g2{x) = Y,hk{b\)Sk, (3.3) 

fc=i 

where the Si are the 6 independent Weyl- invariant scalars ( |B.1| ). The functions^ hi and A 
win be determined in due course by solving the condition 6i = 0. Once this is done, the 
expression for the position of the event horizon will take the form 

rH = lU + J2'^kSk] ■ (3.4) 

The approach developed in [l|, |2| starts by determining a vector field, denoted by v, which 
is normal to the outer marginally trapped surfaces which foliate the horizon A. Two 
properties of v which are essential to its construction are the fact that v is tangent to A, 
and that it must be surface forming, that is, it must satisfy the Frobenius condition 

vAdv = 0. (3.5) 

It will be shown in the following that for the geometry under consideration, up to second 
order in the gradient expansion, these conditions together with Weyl invariance determine 
V uniquely. 

To determine the vector field v one first needs the normal covector to a surface of the form 



(|3.1| ). This is given by m = dS, which up to second order in the gradient expansion is 

m = r db + b dr + dX. (3-6) 

It is convenient to write the normal in terms of the Weyl-covariant derivatives (given in 
(|Aj), (1X1)). One then has^ 



m 



[rV^b + A' {b^V^q + ?,qb^V^b)) dx^" + b{dr + rA^dx^") . (3.7) 



The vector v is to be tangent to A so it must satisfy v ■ m = 0. To solve this condition 
it is convenient to use a special coordinate choice. As discussed in detail in |g], given a 



*In principle these functions also depend on the other Weyl-invariant scalar br, but these formulae are 
eventually evaluated on hypersurfaces r — 1/6 up to terms of second order in gradients. For this reason one 
can set &r = 1 in these functions, since they always appear multiplied by second order scalars. 

^Here A' = X'{b''q). 



parametrization y^ of the horizon so that S{r{y),x{y)) = const, one can choose the gauge 
y = X and then the vector field v takes the form 

_ f}j d _ fdSyW dS\ d\ 

Requiring that the vector v be Weyl invariant imphes that (up to second order) the "bound- 
ary" components v^ take the form 

+ b(aoVo^' + ail^ + Y.'"^^k] +^''I]efc5feJ (3.9) 

\ fc=i / fc=i / 

where a^, c^, e^ are functions of the Weyl-invariant combination b^q. Using (|3.8| ) one finds 
that the r component of v is 

v'' = -brA^uf" - ru^'V^b - br {bA^ + V^b) (ao^o'' + all^') . (3.10) 

It is computationally convenient to normalize v so that 

m2+t>2 = 0, (3.11) 



which implies in particular that coefficients of the longitudinal terms in ( |3.9| ) vanish: e^ = 
for A; = 0, . . . , 6. 

The remaining coefficient functions appearing in v are also not arbitrary. As discussed 
earlier, to ensure that the vector v defines a foliation one has to impose the Frobenius con- 
dition (3.5). As in [Q] the vanishing of v^^dyVp^ is automatic, but the conditions v^^d^Vp^ = 
are nontrivial. In fact these conditions determine the remaining freedom in u. At first order 

in the gradient expansion one finds 

366^2 (56^2 ^ 2) 



oo 



4 (612^4 _ 56^2 _ 2) ' 






while at second order one obtains 



1 



1 66l2fl4^2 



C2 



2 {l + b^q^f 



2V3b^Kq^F,{b\l) 

"' = 6VT1 ' ^ ^' ^ ' 

'^ - (56g2 _ 2)2 + 2(6V-2)(6V + 1) -^^(^'^'^)' 

C5 = -W^5 (ft''?, 1) . (3.13) 



To show that the Frobenius condition (|3.5D is really satisfied one also has to use the relations 

V^Vq,-V,Vq^ = 0. (3.14) 

(valid up to second order in gradients). These relations (discussed further in appendix ^ 
follow from the equations of hydrodynamics. 

This way one finds (as in the uncharged case Q) that the foliation vector v is completely 
determined once A is fixed. As discussed in 0], it seems plausible that this will also be the 
case at higher orders in the gradient expansion. Since the Frobenius condition was imposed 
for the full spacetime (rather than just on the horizon), this vector v actually gives rise to 
foliation of the full spacetime, at least in a neighborhood of the horizon. 

Dynamical quasilocal horizons are spacelike and so m should be timelike and v spacelike. 
Without loss of generality one can assume that m is future pointing and v is outward 
pointing. Then the null normals to the leaves of the foliation of A can be expressed as 

y = £-Cn 
m = i + Cn, (3.15) 

where the scalar C is called the evolution parameter [^, |25| : 

C = \v'. (3.16) 

The sign of the evolution parameter indicates whether A is spacelike (C > 0), timelike 
(C < 0) or null (C = 0). The signs of the coefficients in ( 3.15] ) have been chosen to ensure 



that both i and n are future-pointing, and i is outward-pointing while n is inward-pointing. 
Given the formulae for m and v it is straightforward to calculate C, i and n explicitly. 

To determine the position of the apparent horizon one needs to calculate the null expansions 

9n = r'^aUb , (3.17) 

where the metric induced on the foliation slices is 

Qab = 9ab + Lnb + hna , (3.18) 



where Qat are the coefficients of the metric (|2.7| ) . Using the results of the previous section 
one finds (up to second order) ^"^ 

6 

e£ = 36Sr + ^ef^5fc (3.19) 

fc=i 

o 6 

^- = -T; + T.^n^^k (3.20) 

fc=i 

where the coefficient functions 9^ and 0n are given in the appendix (Q). Note that the 
results are manifestly Weyl-invariant. There is no correction at first order, as required by 
Weyl invariance. With these results in hand, it is straightforward to determine the location 



of the apparent horizon by solving 6(^{rAH) = 0. One finds the result (3^) with 
AH _Ki{b^q,l) 






,AH 



9 - 2bS^ 2612^2^4 



2 19/'/,6„2 _ o\ ^ c /!,fi„9 , i\2' 



h^ 



12 (66g2 _ 2) 5 (56^2 + 1) 

1 
12 (66g2 _ 2) ' 



^^ _ K4 {b^q, 1) b^q'^ {b^^q^ - l^b^^q^ - 2Ub^^q'^ - 232b^q^ - 64) 
4 - 6V-2 ^ 32 (66g2 _ 2)4 (66^2 + i)2 



,^// 



ETslb^'g,!) 6V(bV + 2) 



6V-2 4 (66^2 _ 2)2 (56^2 + 1)' 



^^^ _ i-fe (fe^g, 1) VSb^Kq^ (36i2g4 + 1456^2 ^ 8^ 



^' 6V-2 + 4(6V-2)'(6V + 1)' ■ ^^'^^^ 

In the uncharged case considered in M, u], only hi differs from the result for the event 



horizon. In the present case, the event horizon determined in |19] differs from the apparent 
horizon in the coefficients hi and /14: 

' ^ 2(66g2_2)4 ^ ^ 



The expression 



&(66g2_2)M3 2 (66^2 _ 2)2 M" ^ ^ 



"This computation is fairly lengthy. 



10 



explicitly shows^^ that the apparent horizon lies within (or coincides with) the event horizon 
in the sense that an ingoing radial null geodesic will cross first the event horizon and only 
then the apparent horizon, since r is an afhne parameter on such geodesies. It is also easy 
to check that the apparent horizon is spacelike or null 



1 (l 6V(6V + 2) ' 

(2-6V) U '+ 2 {Wq2-2) 



CirAH) = ,, ,,„,, I ^S, + ^A_i__Z_5, > , (3.25) 



as required. 



4. Hydrodynamic entropy currents 

In hydrodynamics the entropy current is a phenomenological notion constructed order-by- 
order in the gradient expansion. The leading term describes the flow of thermodynamic 
entropy in the perfect fluid approximation. In the case considered in this paper the most 
general form of the entropy current consistent with conformal symmetry up to second order 
in gradients reads 



S" ^ 



4G^ 



+ b Lv,^+hi^ +J2jkvA + (yz^i^^) «") • (4-1) 



The overall factor of I/AGn in (O) comes from the holographic formula for thermodynamic 
entropy. The coefficients j^ and jjl appearing in this expression are functions of the Weyl- 
invariant variable b^q and should be such as to ensure that the divergence of this current 
is non- negative on all solutions of the equations of hydrodynamics. For the case without 



charge this condition was analyzed by [26|, where it was found that it is not possible to 
determine all the coefficients without some additional input. Thus, if the notion of local 
entropy production in the near-equilibrium regime makes sense, there must be some further 
constraints on the form of the hydrodynamic entropy current. This prompted the analysis 
of lH], which showed that the ambiguity in the entropy current is reflected on the gravity 
side precisely as an ambiguity in the choice of horizon used to define the entropy current. 
In [|[] a formula was proposed, which associates an entropy current with each hypersurface 
which satisfies the Hawking are increase theorem, such as the event or apparent horizons: 

3^ = 77^1^-^"^ (4-2) 



"^Both Si and S4 are manifestly positive. 



11 



where the prefactor involving Gm has been introduced to reproduce thermodynamic entropy 
at leading order. This formula, as discussed at length in [|^, is to be evaluated on the chosen 
horizon. The divergence of this current was shown to be proportional to the quantity 
9 1 — C9n, which is non-negative on the basis of the area theorem. 



To apply ([4.2| ) to the gravitational setup of section y one needs the form of v and a 
computation of the determinant of the bulk metric G. As discussed earlier, the vector v is 
completely fixed by the self-consistency of the bulk construction. The right hand side of 
( |4.2| ) evaluated on a hypersurface of the form (3^) leads to ( |4. 1|) with coefficients of the 
first order terms given by 

3^6g2 (56^2 ^ 2) 
Jo - 



4(66g2_2)(66g2 + i)^ 






Coefficients of the transverse second order terms are given by 

1 



1 66^2^2^^ 



32 



2 (56g2 + 1 



,2' 



.3- ^ ^-^^;:f:r'^^ -^3(53.1), 

._ 36V 35V(6V + 2)F,(b3g,l) 

^' ~ (66g2 _ 2)2 + 2 {hS^ - 2) {hS^ + 1) ~^' ^^ '^' ^> ' 

is^ = -^^5 (b^^tZ, 1) • (4.4) 



These are all fixed independently of the hypersurface on which the formula ( [4 .21) is evalu- 
ated. The coefficient functions of the longitudinal terms, proportional to the 4-velocity n^, 
do however depend on the choice of horizon through the functions hj.: 

jI =-F2(6Vl)' + ^^i(fc'g,l)+3/ii, 



2 

612 (|5 _ 12^2) g4 + 1066^2 + 5 
10(66^2 + 1)2 



^2 = ,^,,« ^ , ,,2 + 3^2, 



J3 = 3/i3, 

jj = -L4(6\,l)+3/i4, 

j| = ^Ls (63(7,1) +3/15, 



12 



i = \lq {b\ 1) + 3/i6 . (4.5) 

Evaluating (5^) on the event horizon found in p9|] leads to the result 



II 9 / c! X 3 , o , 3Ki (b^q, l) 1 



2 '^ ^' ^ 6V-2 (6V-2) 



■^2 46V -8' 

II 1 



Js 



466g; 



2 



II _ 3 .3 3K4 (fc^g, 1) 36i2g4 (17612^4 ^ 2866g2 + 20) 

J4-2^4(M,1)+ ^6,2 _ 2 + 32(6V-2)^(6V + 1)' ' 

.11 _ 3feV(bV + 2) ^ 3^ ,^ , 3i^5(^^g,l) 

^^ 4(6V-2)'(6V + 1) 2^^'^ ^^ ^ ^V-2 ' 

.11 3^/369^.^3 (361V + 146V + 8) 3.0 . , 3i^6(6Vl) ,, ^, 

^^= 4(6V-2)^(6V + 1)^ +2^a(M,l)+ ,^,2 _ 2 " (^-S) 



For the apparent horizon one finds the same results as above apart from the coefficient 
functions jf and Jl : 

-UEH) _ MAH) ^ 1 

•^1 -^1 (6V-2f' 

6V(&V + 2)' 



,\\{EH) _ .\\{AH) _ 

2 (66g2_2) 



jr"" - fr"" = 3^ L% ,4 • (4-7) 



In the limit of vanishing charge one can check that the results of [^, ^ are reproduced^^. 

The second order contribution to the entropy current cannot be expressed uniquely in 
terms of the transport coefficients. At first order however there is no ambiguity, and 
indeed the result obtained above is unique at that order. One can easily see that this 
result is consistent with the purely hydrodynamic analysis of Son and Surowka |^, who 
expressed the entropy current at first order in terms of the transport coefficients. This is 
discussed further in appendix ^. 

Calculating the area theorem on either one of the horizons up to the second order in 
gradients one obtains 

uf — Ounl, „ = -. TT. — :tt i^l + 3 ttOa , (4.8) 

nW-VAu (2-6V)l 2 (6V-2)^ / 



^Note however that in the pubHshed version of Q there is a sign error in eq. (81). 
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which matches the hydrodynamic result up to second order in gradients. Thus the di- 
vergence of both entropy currents considered here is non-negative as long as the charge 
remains sub-extremal (i.e. if'q^ < 2). 

Calculating the third order contribution in the bulk requires determining the geometry 
to third order in gradients, which has so far not been obtained. Nevertheless the match 



is guaranteed by the formula (4.2) relating divergence of the entropy current to the area 



theorem on the apparent horizon modulo modifications of bulk-boundary map. 

5. Summary 

The ambiguity in the definition of the entropy current in relativistic second order hydro- 
dynamics was recently connected with the issue of black hole boundaries in the case of 
non-static black objects |l|. This was based on the holographic representation of strongly 
coupled super symmetric Yang-Mills plasma in the framework of fluid-gravity duality Q. 
Reference Q proposed an explicit formula which associates an entropy current with a hori- 
zon in the bulk geometry. As long as the horizon satisfies the area theorem, the divergence 
of the hydrodynamic entropy current is non-negative. In references [Q, ^ entropy currents 
defined by the event horizon and the Weyl-invariant apparent horizon we considered in the 
case of conformal hydrodynamics with no conserved charges beyond the energy-momentum 
tensor. 

This paper continues this line of research by exploring the holographic representation of 
entropy currents in the case of hydrodynamics with a conserved charge. As in [13, , currents 
associated with two horizons were considered here: the event horizon (whose location in the 
relevant geometry was established in |I^) and the Weyl-invariant apparent horizon, which 
was found in this paper using a novel approach to locating apparent horizons developed 
in . This method is particularly suited to situations, where an apparent horizon which 
respects a given symmetry is sought. The work reported here provides another (much more 
complex) example where this approach can successfully be applied. 

The key idea behind Q is that the apparent horizons of interest in the context of fluid- 
gravity duality are only those which are covariant in the sense of the dual hydrodynamic 
description, i.e. can be covariantly specified (in the boundary sense) in terms of b, u^, 
h^jj and their gradients. An important further constraint in the present context is Weyl 
invar iance. 
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One of the main results of this paper is that up to second order in gradients a unique Weyl- 
invariant apparent horizon exists which is covariant in the hydrodynamic sense. As in the 
case without a conserved charge, this horizon is isolated at leading and first subleading 
orders of the gradient expansion and becomes spatial once second order gradient contribu- 
tions are included. As in [|l|, the apparent horizon gives rise to a notion of hydrodynamic 
entropy current which satisfies all the hydrodynamic constraints^^. 

At first order in gradients the entropy current in hydrodynamics with charge was considered 
by Son and Surowka Q. They carefully analyzed hydrodynamic constraints, allowing for 
the possibility of a U{1) anomaly, and succeeded in expressing the entropy current at first 
order in terms of the transport coefficients. As shown in section Q, the result of this purely 
hydrodynamic analysis is reproduced by the holographic entropy current formula proposed 
inQ. 

At second order, in the case of conformal hydrodynamics without charge, a purely hydro- 
dynamic analysis of the allowed form of entropy currents was performed in ||2^ (see also 
p!^). The current associated with the apparent horizon was found to be consistent with 
that analysis |j^]. In the case with charge such an analysis has not yet been carried out; it 
would furnish an nontrivial check on the results obtained in this paper. 

It would be very interesting to study apparent horizons in the cases of non-conformal |27|, 



and super fluid |28, 29 1 fluid-gravity dualities and explore the issue of entropy currents 
in those contexts. Another possible direction of further research would be to include 



background flelds along the lines of ^^, |30| 
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A. Weyl covariance 

Conformal symmetry of AA = 4 super symmetric Yang-Mills theory can be extended to the 
bulk as follows [Isl, ^: 



gi^u -^ e-'^^g,,^, u^ -^ e'^u^, h -^ e'H and r -^ e'^r (A.l) 

where ^ depends on the coordinates x^ [23]. A quantity which transforms homogeneously 
with a factor of e""'^ is said to transform with Weyl weight w. 

A beautiful formalism allowing for manifest Weyl covariance in conformal hydrodynamics 



was introduced by Loganayagam [11 1 and applied to fluid- gravity duality in [p^ , p3| . The 
basic tool is the Weyl-covariant derivative 2?^, which preserves the Weyl weight of the 
differentiated tensor. It is constructed using the vector field Av defined by [|ll| 



Au = u^VxUu )—u^ . (A.2) 

O 

This quantity is of order one in the gradient expansion and transforms as a connection 
under Weyl-transformations 

Ay^ Au + d^cf) . (A.3) 

Due to this property it can be used to compensate for derivatives of the Weyl factor when 
differentiating a Weyl-covariant tensor. For instance, one has 

V^b = d^b - A^b, 

'^t.q = df,q + 3Af,q . (A.4) 

For further details the reader is referred to the original literature cited above. 
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B. Lorentz tensors at second order 



Scalars: 



Tensors: 



Si = b^a^.a^'', 

S3 = b^n, 

54 = b^-^P^^'^V^qV.q, 

55 = b'^q-^P'^^V^V.q, 
Se = b''q-^Pf'''LV,q. 



fB.l) 



Vectors: 



Vi. = bP.uV.a'^P, 



V2^ = bP^^VpCO^P, 



V, 



3fi 



bl^a^x, 



V4f, = bq ^a^'^Vaq, 



1, . ctn 



V5^ = bq''u^''V^q. 



(B.2) 



■ IflU 



''2fiu 



^2, 

Tsuu 



uf"Dpa^u, 



1 



a/3 



Ti 



lO^j/ 



Til 



n^^^a^/3, 

J^ a/3 ^ X cr 

2^ X{fi^al3iy)aU U , 

q-^Uf^V^qV^q, 

q~^Uf^V^V^q, 

q-'ufj^Vpq, 

al3X 



^ apX -l-n 

iMi^ = 2^(m ^'^)^^"9 ^/39- 



(B.3) 
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Here li'^u is the projector which can be used to create symmetric, traceless tensors: 






Up^p^^+p^pf^-^p-^p, 



flU 



The scalar TZ is defined as in ||2^ and C^aui3 denotes the Weyl tensor: 



C, 



fiuXa 



R 



1 



puXa 



d-2 



[9n[\Ra]v - 9u[xRa]fj) + 



1 



(d-l)(d-2) 



9ti[x9a]iyR ■ 



(B.4) 



(B.5) 



C. Second order hydrodynamics 



As shown in [O] the vahdity of the gradient expanded solution requires the conservation 
of 



1 



1 + 6^ 



(PuU + ^UnUu) 



2n^, 46V(-l + ft'g'(12Ac2-l)) 



2(7^1/ 2(1 + cijTi^i, 2T2^i, 2ci T^^^y 



+ 



62 



+ 



62 
2^/Zh'^q^K 



1 + 



^5m!^ + ]^ _^ ^6^2 -^6^(1^ + ^2'^8ij.u 



62 



62 






(C.l) 



and 



J„ 



1 fV^qUf, Sb^q'^Klf, V^b^q{2 + b'^q'^) 



^ttGn V 2 



+ 



+ 



+ 



3\/36'^g^K2 



^2, 



2(l + 66g2; 
364^^2 



:i + 66^2)^2 



Vc 



O^t 



3\/36g 
(1 + 66g2) 



Via + 



lY'^lM 



2a4(6V + l)+V36V,, , 



vl + 6V)2'"'' 2(66(72 + 1)2'^^ ' 1662(66^2 + 1)2 

as (66g2 + 1) + ^^9 (24^2 - l) q^ - ,/3b^q 



862 (^6^2 + 1)^ 



-V5^, 



(C.2) 



where all the quantities appearing above are given in |19]. 



This leads to the equations of hydrodynamics 



V„b 



67gVo^ , bVi^ , 367gV4^ , Siu^ 



+ 



+ 



+ 



b^'q^SjU^ 



2-6V 6V-2 (6V-2)^ 6-366g2 4 (66^2 _ 2) (66^2 + i 
66^2 (56^2 + 2) S5U^ 2V3b^Kq^Seu^ 



+ 



4(66g2_2)(6V + l) (6V - 2) (6V + 1) ' 
^^ , b^q^S^u^ q{b\^ + 2)S^u^ 2V^b^Kq^S^u^ 



466 g2 + 4 
From these it follows that 



Onj\.it OpJ\ 



^fl^U ^U^fl 



4 (67g2 + b) 
b-'q^ 



b^q^ + 1 



-XVyVQ^-V^V^y), 



(C.3) 



-66g' 



(C.4) 



which leads to the relations (3.14) used earlier. 
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D. The expansions 

The coefficients appearing in the expansion Op are^^: 
9} = SbrKi - 2b^Br^F2F2' + -b'^Br^Li', 

' ~ 465^7 + 5^9 (56^2 + 1)2 ^^4^2 _ 52^4 _ ^2) I 

+ lOfe^OgV^O - 56i8^V8 - 5616^^6 + 5b''q\^ {q^ + 2r6) + 5b'^q\^ {q^ - Ar'') + 

- 46^° (2g^ + 5gV^) + 126^^V^ + 96^ {2q^r'^ - r^) - b* {Sq^ + Or^) + 136V^ + 13r 
1 



^' "46^' 



^4^ »V / 

' 32r (&6q2 _ 2)3 (56^2 + 1)2 (1 _ 52r2)(_56^2 + ^2^2 + 54^4) V 

mW^ {b%^ - 2) {b\^ + 1) {b\^ + 2f- 126r3 (66^2 ^ 1^2 ^^ig^g ^ 14^12^4 ^ 12^6^2 ^ g^ ^ 

- 326^6 (56^2 ^ 1^2 ^^12^4 _ 7^6^2 _ 2) + i6fe4 (56^3 ^ ^^2 ^^12^4 ^ 8^6^2 ^ 4^ ^ 

- r2 (56^2 (9524^8 ^ 2146^^9^ + 5806^2^4 ^ 45556^2 ^ g4^ ^ g4^ ^ 
+ 126V^ (66g2 ^ 1^ (^^18^6 ^ 14512^4 ^ 1256^2 ^ g^ ^ 

- 125^ g2^ (^56^2 ^ 1^ (^^18^6 _ 18512^4 _ 2056^2 + 8) ) + 

4 (66 g2 _ 2) (66^2 + 1) (1 _ 52^2) (_56 ^2 + ^2^2 + ^4^4) \" ^ O q ^ 

+ 26^6 (g6 _ 5/^6^ ^ g^lO (^4 _ ^2^6) _ g^6^2^2 ^ 12^5^2^ ^ 4^4 (^2 ^ ^6) _ 4^2^ ^ 
, ^M^V±llj^ , 3 2 , 6V(fc6g2_3fc4^4^i)^^2 

+ c,orA4+ 56^2 _ 2 a(63g) + 2 ^ ^ ^ 2(1 - 62r2) (-66^2 + 52^2 + 54^4) ' 

.1 3 n n , 65g2 (3^7 2^ 4^6 2 _^ g^^ _^ 4\ 

^f = U^l>^ f 2fe23^6^7 ^17 4^ (2 g^6) 4^16 4^6 _ 3^15 4^5 

' 4r6 (66^2 _ 2) (56^2 + 1)2 (^4^2 _ 52^4 _ ^2) V 'i \^ ) H 

_ 46l4(?2rl0 - 3613<;V3 - 4612^2^8 ^ 4^11^2^7 ^ 4^10^2 (^2 ^ ^6) _ g^9^2^5 _ 4^8^4 (^2 ^ ^6^ ^ 

- 66^g2^3 _ 4^6^2 (^^2 _^ ^6^ _ 4^5 ^2^ _^ 4^4^2 _ 4^2^4 _ 4^2\ _^ 3^^^^ _^ -b'^Br'^Le' . (D.l) 
^*The prime denotes a derivative with respect to br. 
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For On the result is 



^i = 2F2F2' - \l,\ 



1 662^2^4 



q2 I 

" 6V3 5^7 (56^2 + 1)2 ^_^4^2 + ^2^4 + ^2)2 V 

6^(96^9^ + mb^'^q^r^\b^q^ + 1) + r^(56^2^^ _ 2856^2 _ g) ^ 
+ b'^r^bbWb^'^ - 2) - 6) + 26^^2(566 g2 _ 4) _ Qh^r\bb\^ + 2)) ) , 

4 _ fe'^gV^ 



'n 



16 (66g2 _ 2)'^ (66g2 + ly (52^2 _ 1)^ (_^,4g2 + ^2^4 + ^2)^ 

ISft^gV^ {b\^ - 2) (6^2 + 1) (56^2 ^ 2)2 + i26r3 (fo^^^ ^ ^^2 ^^ig^g ^ ^4^12^4 ^ ^266^2 ^ g^ ^ 
+ 166^*5 (6^2 + 1)2 (512^4 ^ 8^6^2 ^ 4^ ^ 3254^2 (^6^2 ^ 1^2 ^2^12^4 ^ ^6^2 ^ 2) + 

- r^ (56^2 (57524^8 ^ 2626^^9^ + 4366^2^^ ^ 2166^^^ - 32) + 64) + 

- 126V^ {b^'q^ + 1) (6^^^ + Ub^^q^ + 126^ g^ + g) + 

- 1265gV (66g2 ^ 1) (3518^6 ^ 10^12^4 ^ 4^6 ^2 + 24) ) + 



^2„4 77 

'^ '^ ^0 '' ^^(.4 2/0 , n^6 2 , nul2i^ 



26V(2 + 96'r + 76^''7^) + 



2 (66^2 _ 2) (66^2 + 1) (56^2^2 _ 54 (^2 + ^6) + ^2)2 

+ 365g2(_4 ^ 5i2^4y _ 254(2 + 356^2 ^ ^12^4)^6 ^ 2(2 + hb'^q^){r + 6^9^)2) + 



53^5 (256^2^2 _ 354^2 _^ 2r2) Fo^ 265^r'^ (56^2 _^ 1) Qp^ 3 ^ 

(56^2^2 _ 54 (^2 + ^6) + ^2)2 + (66^2 _ 2) (62r2 - 1) (64^2 _ ^2^4 _ ^2) ^(^58) " 2^^ ' 

5 3 , Fo b^q^ {-3b'^q'^r + 4fe^g2 _ g^^ _^ 4) 



2 ^ 262Br2 8Sr3 (6V - 2) (6V + 1) 



^55 _3 



""^ ^66i5(z'r^ + 3b'\\^ + 46i2<^V2 + Sft^V^ + 



2r2 (66g2 _ 2) (66^2 + 1) (_54^2 + 52^4 + ^2)2 
+ Ub^q'^r^ - Ub^q'^r* + eft^gV^ + b^ (l2r^ - SgV^) - 46^ (g^ _ 6r'') + 246V^ + I2r'^) + 

Kg 



+ 



^56 _3 



2r2 (66g2 + 1)2 (52^2 _ 1) (_54^2 + 52^4 + ^2)2 

12r (-ftlOg^ ^ 58^2^4 ^ 56^2^2 _ 54^2 ^ 52^4 ^ ^2^ ^^ ^ 
+ 63g2 (^54^4 ^ 52^2 ^ 1^ (357^2^ _ 456^2 + g5^ _ 4) \ _ 3^g/ ^ p_2) 
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E. The entropy current at first order 

The authors of reference [^ considered the most general form of entropy current and 
constitutive relations in hydrodynamics with a charge current (allowing for the possibility 
of a C/(l) anomaly). They showed that at first order in gradients the entropy current must 
be of the form 

S^' = su^" - ^v^' + -Dl^" (E.l) 

where s is the equilibrium entropy density, D is a coefficient discussed below, and v^ is 
first order correction to the charge current: 

J/^ = n-u^ + u^" (E.2) 

By symmetry arguments this correction can be expressed as 

z./^ = -^aP^-a. (I) + ie/^ (E.3) 

where a and ^ are transport coefficients. For the fluid considered in this paper these 
transport coefficients can be read off from (C.2). For this one needs the relations 



r.^. „.V1^ ,E.4) 



from which one derives 



vr; f6V-2)^ ^ ' 



The coefficient ^ is non- vanishing only if the theory is anomalous (in the sense that coupling 
it to a background gauge fields results in an anomalous divergence of the charge current). 
As shown by Son and Surowka Q , ^ is proportional to the anomaly coefficient C: 

This equation allows one to express C in terms of known quantities, and then one can 
compute D using the relation ^ 

D = -C^ . (E.7) 

ST ^ ' 



Proceeding in this fashion one can express the entropy current ( |E.lD in terms of the variables 
b and q used in this paper. This leads to 



1 /l ^ 36V(6V + 2) ^^, Vsb-^ 



Kq 



3 



AGN\b3 4(66g2_2)(6V + l) ° (6V + 1) / 

This is precisely the first order part of (j^^) with coefficients ( [4. 3D , which were derived from 
the holographic entropy current formula ([4. 21) proposed in ^. 
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